We present methodology to derive high-precision estimates 
of the fundamental parameters of double-lined spectroscopic 
binaries. We apply the methods to the case study of the double- 
lined (3 Cephei star (3 Centauri. We also present a detailed anal- 
ysis of /3Centauri's line-profile variations caused by its oscil- 
lations. High-resolution spectral time series and visual or inter- 
ferometric data with a good phase distribution along the or- 
bital period are required. We point out that a systematic er- 
ror in the orbital amplitudes, and any quantities derived from 
them, occurs if the radial velocities of blended component lines 
are computed without spectral disentangling. This technique 
is an essential ingredient in the derivation of the physical pa- 
rameters if the goal is to obtain a precision of only a few 
percent. We have devised iteration schemes to obtain the or- 
bital elements for systems whose lines are blended through- 
out the orbital cycle. We derive the component masses and 
dynamical parallax of (3 Centauri with a precision of 6% and 
4%, respectively. Modelling allowed us to refine the mass es- 
b— h ' timates to 1% precision resulting in M\ — 10.7 ± 0.1M Q 
^5 , and Mi = 10.3 ±0.1 M , and to derive the age of the sys- 
' tem as being (14.1 ± 0.6) x 10 6 years. We deduce two oscil- 
lation frequencies for the broad-lined primary of /3 Centauri: 
/i = 7.415cd _1 and /2 = 4.542cd _1 or one of their aliases. 
The degrees of these oscillation modes are higher than 2 for 
both frequencies, irrespective of the alias problem. No evidence 
of oscillations in the narrow-lined secondary was found. We 
propose that our iteration schemes be used in any future deriva- 
tions of the spectroscopic orbital parameters of double-lined 
binaries with blended component lines to which disentangling 
can be successfully applied. The combination of parameters 
resulting from the iteration schemes with high-precision esti- 
mates of the orbital inclination and the angular semi-major axis 
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solution of the system. 
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1. Introduction 

Binary stars have long been considered as astrophysical labora- 
tories, providing one of the best tests of stellar structure models 
(e.g., Maceroni 2005 and Ribas 2005, for recent reviews). This 
is particularly so for massive binaries because their structure 
and evolution are not well-understood, while being of great im- 
portance for the chemical enrichment and evolution of galax- 
ies. The convective and rotational mixing properties of mas- 
sive stars with a well-developed convective core are still poorly 
calibrated, while being the dominant factors determining their 
evolution (e.g. Maeder & Meynet 2000). While observational 
capabilities and analysis tools have improved significantly in 
recent years (e.g. Hilditch 2004a), there is still a lack of ultra- 
precise fundamental parameter determinations of binaries with 
an OB-type component (e.g. Hilditch 2004b). Indeed, compo- 
nent mass estimates with a precision better than 2% are avail- 
able for relatively few such systems, although such a preci- 
sion is necessary to provide stringent observational tests for 
stellar structure and evolution models (e.g. Andersen 1991). 
In this paper, we provide methodology to achieve a high pre- 
cision for mass estimates from combined interferometric and 
spectroscopic data of double-lined spectroscopic binaries with 
merged component lines, and we apply it to the massive bi- 
nary /3Cen. The methods are based on spectral disentangling 
(Hadrava 1995, 1997, 2001, 2004b). 

The bright star [3 Cen (HD 1 2245 1 , HR 5267, B 1 III, m v = 
0.6) has been known to be variable in velocity since the be- 
ginning of the twentieth century. It is the brightest member of 
the class of (3 Cephei stars, a homogeneous group of oscillat- 

* Radial-velocity data available electronically from the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) 



ing B0-B3 stars (see Stankov & Handler 2005 for a recent re- 
view). They have low-degree, low-order pressure and gravity 
modes with periods of a few hours excited by the k mechanism 
(Pamyatnykh 1999). They reveal amplitudes of several tens 
of mmag down to the detection threshold in UBV, so these 
stars are good potential targets for in-depth seismic studies. 
Asteroseismology of j3 Cephei stars indeed received a lot of 
attention lately, after it became clear that their oscillations can- 
not be explained in terms of standard evolution models. For 
two prototypical class members, the oscillations revealed dif- 
ferential internal rotation and the occurrence of core convec- 
tive overshooting (HD 129929: Aerts et al. 2003 and Dupret et 
al. 2004; i^Eri: Pamyatnykh et al. 2004 and Ausseloos et al. 
2004). 

As one of the brightest stars in the southern hemisphere as 
a whole, (3 Cen has been the subject of numerous studies. We 
refer to Ausseloos et al. (2002, hereafter Paper I), Davis et al. 
(2005, hereafter Paper II), and references in these two papers 
for an overview of these studies, without repeating all of the 
information here. We summarise only briefly the characteristics 
of the system that are relevant for our current work. 

High-resolution spectra covering 12 years revealed that 
/3Cen is a double-lined spectroscopic binary with an orbital 
period of 357 days and an eccentricity of about 0.81 (Paper I). 
Interferometric data assembled with the Sydney University 
Stellar Interferometer (SUSI) and covering 7 years led to sim- 
ilar values for the period and eccentricity and, moreover, to 
an orbital inclination of 67.4° and an angular semi-major axis 
of 0.0253", as well as to a brightness ratio of 0.868±0.015 
(Paper II). The spectroscopic variability is due not only to the 
binarity, but also to oscillations of the components with periods 
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of several hours (Paper I). The star is photometrically constant 
at the level of mmag. 

The independent orbital fits to the spectroscopic and inter- 
ferometric data have four parameters in common: the orbital 
period (P rb), the epoch of periastron passage (-Bo), the ec- 
centricity (e), and the longitude of periastron (u>). The values 
of these four differed by less than the ler uncertainties such 
that a fully consistent orbital solution for (3 Cen was achieved 
(Paper II). The combination of all the available information 
subsequently led to the conclusion that the system consists of 
components with equal masses of Mi = M2 = 9.3 ± 0.3 M@ 
(i.e., a precision of 3.2%) and that it has a dynamical parallax 
of 9.78±0.16mas. 

We show in this work that the masses of the components 
of (3 Cen were significantly underestimated due to a system- 
atic error in the amplitudes of the spectroscopic orbit. This is 
a consequence of the inappropriate way in which the radial ve- 
locities were estimated from line profiles of merged spectra. A 
similar conclusion was reached recently by Tango et al. (2006) 
for the triple system A Sco and occurs for any spectroscopic bi- 
nary in which both components contribute to the lines used for 
the orbital radial-velocity (hereafter abbreviated as RV) deter- 
mination. Earlier attempts to avoid such systematic error can 
be found in Tomkin et al. (1995) for the i5 Scuti star 9 2 Tauri. 
In that work, the authors subtracted the lines of the primary 
by means of spectra of reference stars with the same spectral 
type before computing the secondary's RV values. We provide 
analysis schemes based on spectral disentangling to overcome 
this problem of systematic errors in a more accurate way. Our 
schemes allow us to eliminate the systematic errors in the phys- 
ical parameters. We illustrate our method by its application to 
the case of (3 Cen. Our methodology is applicable to the analy- 
sis of any spectroscopic binary whose line profiles can be suc- 
cessfully disentangled. It leads to a significant improvement in 
the precision of the physical parameters and dynamical parallax 
of such systems, similar to the case of binaries with emission- 
line stars (Harmanec 2002). 

2. Methodology for orbital determination 

The data used to illustrate our methodology are the Silll 
A 4552.6 A line profiles of /3Cen obtained over 12 years with 
the ESO CAT telescope and with the Swiss Euler telescope, as 
described in Paper I. By means of illustration of the system's 
lines and of the occurrence of a systematic error in the orbital 
RV determinations used in Papers I and II, we recall here in 
quite some detail the way the RVs were obtained. We show 
in Fig.^a typical line profile of the system averaged over one 
night of data, i.e., a profile in which the oscillatory variations 
are averaged out. As can be seen, the lines produced by the 
two components are blended with each other. In none of the 
available spectra are the two components' lines well separated. 
Moreover, the oscillations induce deviations from a Gaussian 
shape for each of the line components (see Fig. 1 in PapeiT). 
For this reason, the RV values of the component with the nar- 
rowest line were derived from the line centroid (first moment, 
see Aerts et al. 1992 for a definition) in Paper I, with the inte- 
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Fig. 1. Silll A 4552.6 A line profile obtained by averaging ten 
spectra measured on 16 May 1988 (solid line). The two vertical 
lines denote the integration limits which were used in PapeiT to 
calculate the RV of the secondary. The two dotted lines denote 
the disentangled line profiles of both components. 
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Fig. 2. Silll A 4552.6 A line profile obtained by averaging 39 
spectra taken during a period of 12 nights (3-14 August 2000). 
The vertical lines demonstrate the estimate of the total line 
width of the broad-lined component. 

gration limits indicated in Fig.[2corresponding to a normalised 
flux value of 0.95. 

The orbital RVs for the broad-lined component could not 
be derived in this way. The following strategy was therefore 
followed in PapeiT. The full width of the line was derived from 
epochs when the radial velocities of both components did not 
differ much (see Fig.|2j. It was assumed that this width is con- 
stant, which is a reasonable approach because rotational broad- 
ening is dominant for this component. The centre of the broad 
line, obtained from averaging spectra over a night (CAT) or 
over two weeks (Euler), determined by starting from either its 
left or its right wing, was taken as a good estimate of the RV of 
the broad-lined component. 

We show below that these procedures lead to an underesti- 
mation of the true RV values, particularly for the narrow-lined 
component. We provide the final RV values of both components 
in Table 1 (only available electronically from CDS). Hereafter, 
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we will refer to the star producing the deeper and narrower 
Silll lines as the secondary and to the other component as the 
primary. Although this is the opposite of what has been done 
in the literature so far, we show that the component with the 
broader lines is indeed the more massive of the two. 



2.1. korel disentangling 

Starting from the orbital solution presented in Paper II, we 
applied KOREL spectral disentangling (Hadrava 1995, 1997, 
2001, 2004b). Although KOREL was not developed to treat line- 
profile variations due to oscillations, Harmanec et al. (2005) 
showed that the code is able to treat such a complex combina- 
tion of variability. 

KOREL was applied to our /3Cen CES spectra for many 
different sets of code input parameters and weights, the 
CORALIE spectra being too noisy to allow convergence. The 
resulting disentangled profiles were evaluated each time by vi- 
sual inspection, paying attention to smoothness, symmetry, and 
the residual sum of squares. We considered both the situations 
where line strengths were allowed for and were not taken into 
account. Moreover, we used several types of weights. Weights 
proportional to (S/N) 2 turned out to lead to the most stable 
solution. The best results were clearly obtained when KOREL 
was allowed to search through a larger subregion in orbital pa- 
rameter space than indicated by the uncertainties obtained in 
Paper II. 

The adopted KOREL solution was obtained in four subse- 
quent steps, in which the solution of a particular step was used 
as the initial guess in the next step. First the light intensities 
were kept fixed, next the intensity of the primary was allowed 
to vary, then the intensity of the secondary could vary, and fi- 
nally the intensities of both components were allowed to vary. 
This led to a final KOREL orbital solution which was slightly, 
but significantly, different from the one in Paper II, the largest 
discrepancy occurring in the value of the semi-amplitude of the 
secondary's orbit K-x: 63.8 ± 0.6 km/s (PapeiTI) versus 68.1 
km/s (KOREL). The present version of KOREL unfortunately 
does not provide errors in the orbital parameter values given as 
output. 

The dotted lines in Fig^show the disentangled profiles of 
the primary and secondary shifted to the orbital RV obtained 
in Paper I. It is obvious from Fig.^that the RV estimate of the 
secondary derived in PapeiT is an underestimation of the true 
RV and that we must take this into account in the derivation of 
the physical parameters of the components. 

In principle, a single application of KOREL disentangling 
to the blended line profiles should be enough to obtain the fi- 
nal orbital solution. In practice, however, tiny changes in the 
KOREL input parameters changed the final orbital solution con- 
siderably while producing only small changes in the mean dis- 
entangled profiles and rms value. For this reason, we devised 
the iteration schemes discussed below. 
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Fig. 3. A Silll A 4552.6 A line observed on 15 March 1998 (full 
line) and the secondary's disentangled line profile, which is 
shifted according to its corresponding orbital RV (dashed line). 
The residual spectrum obtained by subtracting the shifted dis- 
entangled line profile from the original spectrum is shown as 
dots. 

2.2. Preliminary update of the orbital solution 

The disentangled line profile of the secondary allows us to im- 
prove the orbital solution of the primary. We carried out the 
following procedure for each of the spectra. We shifted the 
secondary's disentangled profile according to its orbital RV ob- 
tained in Paper I. We subsequently subtracted this shifted dis- 
entangled profile from the original blended line profile. This 
procedure is illustrated in Fig. [3] for one measurement. In this 
way, we obtain a A 4552.6 A line profile of the primary whose 
position and shape is due to the orbital velocity, as well as to 
the variability due to the oscillations (the dotted line in Fig.[3j. 
This procedure resulted in 402 line profiles (A, /(A)) which 
were used to calculate the true RV of the primary. In view of 
the bumpy profiles, we cannot use a Gaussian fit to compute 
this true RV. We determine it as follows: 



Wrad = 



/(l-J(A)HdA 



EW 



with EW denoting the equivalent width of the line profile 



EW 



(l-7(A))dA, 



(1) 



(2) 



and v\ derived by the equation 

A -A 

v\ = — : c, 

Ao 

with c and Ao the speed of light and the laboratory wavelength, 
respectively. What we call the true RV is thus the centroid of 
the line profile (also termed the first moment, see Aerts et al. 
1992 for further extensive discussion of this quantity). Its com- 
putation was done by fixing the integration limits interactively 
for each separate profile, after subtracting the disentangled pro- 
file of the secondary. In this way, 402 RV values of the pri- 
mary were obtained (compared with only 27 values used in 
PapeiT obtained with the procedure described above by means 
of Fig. 2). 
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Fig. 4. The disentangled Silll A 4552.6 A line profiles of the 
primary (full line) and secondary (dashed line). 



Next, these 402 RV values of the primary were added to the 
402 values for the secondary derived in PapeiT, and the code 
FOTEL (Hadrava 1990, 2004a) was applied to this combined 
RV dataset. The results are quite similar to the ones listed in 
Paper II, except that the semi-amplitude of the primary's orbit, 
K\, is significantly smaller. This is consistent with the KOREL 
disentangling analysis, which led to a higher value of the semi- 
amplitude of the secondary, K2 (compared to the value found 
in Paper II). Both these results, i.e., a lower K\ value from KO- 
REL and a higher K2 value after disentangling the secondary's 
profile, suggest that the mass ratio M1/M2 has been underesti- 
mated in Paper II. 

Subsequently, a second attempt was undertaken to dis- 
entangle the spectra with KOREL by searching through a 
neighbouring subregion of the orbital parameter space cen- 
tred around the updated orbital solution with a higher mass 
ratio. Again, satisfactory results were obtained similar to the 
orbital solution from FOTEL, except for the values of the semi- 
amplitudes K\ and K^- In comparison with our first attempt to 
apply KOREL, the residual sum of squares was lower. Figure^ 
shows the best disentangled Silll A 4552.6 A line profiles of 
both components at this stage of the process. 

2.3. Analysis of the systematic error 

We return to Fig.^ The merged line profile, /(A), is obviously 
the sum of the primary's average line profile ii(A) and the 
secondary's average line profile /2(A). The primary's disentan- 
gled line profile delivers a good approximation of I\, while the 
one of the secondary results in an approximation of I2 (dot- 
ted lines in Fig.Q. The 402 RV values for the secondary de- 
rived in Paper I were obtained by calculating the centroid of 
the composite Si III A 4552.6 A line profile, i.e., from Eq. Q 
with — ii(A) + /2(A). To minimize the interfering influ- 
ence of the primary's line profile 2j (A), the integration limits 
in Eqs. Q and (|2ji were fixed corresponding to a flux value 
of 0.95 in PapeiT, as indicated in Fig.^ The secondary's 
RV is, however, given by the centroid of only the secondary's 
line profile /2(A). So, calculating the secondary's RVs as was 



done in PapeiT, gives each v\ in the integral Q too high a 
weight. When both components have a similar RV, the method 
adopted in Paper I provides a reasonable approximation of the 
secondary's RV for the following reasons: 

1 . the depth of the secondary's line profile has a significantly 
larger contribution to the depth of the composite line profile 
than the depth of the primary's line profile; 

2. the depth of the primary's line profile does not vary much 
within the integration limits because of the large width of 
the primary's line profile on the one hand and the fact that 
both components have a similar radial velocity on the other 
hand; 

3. the equivalent width in the denominator of Eq. Q nor- 
malises the weighted integral in the numerator. 

A systematic error is, however, introduced when the com- 
ponents' RVs are significantly different. This is the case in 
Fig-H as the spectrum shown is observed close to a time at 
which the primary (secondary) has its minimum (maximum) 
orbital RV. Indeed, in such cases, the second condition above is 
not fulfilled, as the depth of the primary's line profile changes 
considerably within the integration limits. One can derive from 
Fig. n that, within the integration limits, the depth of the pri- 
mary's line profile decreases monotonically with increasing 
wavelength so that the additional weight given to each v\ in 
Eq. Q strongly varies due to the blending of both line profiles. 
As lower v\ values systematically get more weight than higher 
v\ values, the procedure applied so far still underestimates the 
RV of the secondary. We cannot but conclude that a systematic 
error has been introduced in PapeiT in the calculation of the 
RVs, due to the fact that the line profiles of both components 
are so strongly blended with each other. 

2.4. Iterative determination of the orbital parameters 

Since the orbital parameters common to the fit to both the spec- 
troscopic and interferometric data were in agreement with each 
other, it is appropriate to assume that this spectroscopically de- 
termined orbital solution is already close to the true orbit and 
can hence be considered as a good initial solution to start an 
iterative process to improve the orbital parameter values. Two 
different iteration schemes were applied. 

2.4.1. Iteration scheme I 

Figure 5 summarises the features of the first iteration scheme. 
In each iteration step, a new orbital solution is calculated by 
applying FOTEL to the combined dataset of new RV values 
of both the primary and secondary. The new set of the pri- 
mary's RVs is obtained by calculating the centroid of the "sec- 
ondary subtracted line profiles". The latter refer to the spectra 
obtained by taking the difference between the original spectra 
and the shifted secondary's disentangled line profile (abbrevi- 
ated as "secondary's LP" in Fig. 4). The shift corresponds to 
the orbital velocity value given by the orbital solution of the 
previous iteration step. This procedure to determine the sec- 
ondary subtracted line profiles is illustrated in Fig.[3] A new 
set of the secondary's RV values is calculated in each iteration 
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Fig. 5. Flowchart diagram of Iteration Scheme I (see text for 
details). 



Fig. 6. Flowchart diagram of Iteration Scheme II (see text for 
details). 



step in an analogous manner by means of the primary's disen- 
tangled line profile (abbreviated as "primary's LP" in Fig. 4). 
The primary's and secondary's disentangled line profiles men- 
tioned above are properly normalised versions of the best dis- 
entangled profiles obtained with KOREL. All centroid velocity 
values were made with fixed integration limits to reduce the 
noise level. The widths of the primary's and secondary's disen- 
tangled line profiles provide us with an objective way to select 
these fixed integration limits. 

2.4.2. Iteration scheme II 

The second iteration scheme is similar to the first one, but 
slightly more complicated. For the sake of clarity, the reader is 
advised to refer to Fig. 6 while reading the following descrip- 
tion. 

In each iteration step, a new orbital solution is again de- 
rived on the basis of the combined dataset that consists of the 
newly derived RV values of both components. Just as in itera- 
tion scheme I, the new primary's (secondary's) RVs are derived 
by calculating the centroid of the "secondary's (primary's) sub- 
tracted line profiles" with fixed integration limits according to 
the width of the disentangled line profile of either component. 



The difference between the two iteration schemes lies in the 
way that these "primary's/secondary's subtracted line profiles" 
are constructed. For each original spectrum, both schemes shift 
the primary's (secondary's) disentangled line profile according 
to the corresponding orbital velocity that was found in the pre- 
vious iteration step and subtract it from the original spectrum. 
While iteration scheme I uses the orbital RV corresponding to 
the orbital solution derived on the basis of the old RV datasets 
of both components, iteration scheme II derives the new pri- 
mary's (secondary's) RV dataset by using the orbital solution 
derived on the basis of the old secondary's (primary's) RVs. 

Figures0and|8]show the results of the application of itera- 
tion schemes I (open symbols) and II (filled symbols). It reveals 
only minor differences in the final parameter values obtained 
with both schemes. Scheme II leads to a slightly higher value 
of the orbital period, but the differences in the system veloc- 
ity as well as the semi-amplitudes of both components' orbit 
are totally negligible. The stability of the semi-amplitudes K\ 
and K2 is of particular importance as these parameters allow 
an accurate mass determination of both components. The up- 
per panel of Fig. [8] shows that both iteration schemes make the 
eccentricity converge to a value in excellent agreement with the 
interferometric value. One can derive from the middle panel of 
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Iteration step 

Fig. 7. Comparison of the results of iteration scheme I (open 
symbols) and II (filled symbols): evolution of the orbital period 
-forb [upper panel], system velocity u 7 [middle panel], and the 
semi-amplitudes K\ (dots) and K2 (triangles) [lower panel]. 
The vertical bars denote the errors provided by FOTEL. 



Fig.[8]that iteration scheme II puts the value of the epoch of pe- 
riastron passage closer to the centre of the interval which was 
derived for this parameter on the basis of interferometric data 
than scheme I. The lower panel of Fig.[8] reveals an increase 
in the value of the longitude of periastron which eliminates 
the small discrepancy which appeared when applying iteration 
scheme I. 

2.4.3. Evaluation of both iteration schemes 

The final orbital solutions derived with iteration schemes I and 
II are listed in the second and third columns of Table|2] respec- 
tively. 

To reveal the origin of the small differences between the 
results obtained with the two iteration schemes, we examined 
the dataset consisting of the differences between the final RV 
values derived with iteration schemes II and I: 

{Av(t)} t = {Urad,Il(*) - «rad, I (*)}*< 

This dataset includes two subsets: the primary's and sec- 
ondary's radial velocity differences between schemes II and I: 



Iteration step 

Fig. 8. Comparison of the results of iteration scheme I (open 
symbols) and II (filled symbols): evolution of the eccentricity e 
[upper panel], epoch of periastron passage E [middle panel], 
and the longitude of periastron ui [lower panel]. The two hori- 
zontal, dotted lines in each panel indicate the corresponding in- 
terval estimated from the interferometric measurements taken 
from Paper II. The vertical bars denote the errors provided by 

FOTEL. 

Table 2. Orbital parameters for (3Cen obtained from KOREL 
disentangling and by the application of iteration processes to 
the observed Silll A 4552.6 A line profiles listed in Paper I. The 
errors are la estimates resulting from FOTEL assuming the KO- 
REL disentangled profiles to be error-free. 



Parameter KOREL 



Scheme I 



Scheme II 



P orb (days) 356.94 
v-y (fans -1 ) - 
Ki (kms -1 ) 57.4 
K 2 (kms -1 ) 72.3 
e 0.825 
E (HJD) 



356.86 ± 0.03 356.92 ± 0.03 

9.1 ±0.3 9.3 ±0.3 

63.3 ±0.6 63.2 ±0.6 

72.7 ±0.7 72.1 ±0.6 

0.824 ± 0.002 0.825 ± 0.002 
2451600.08 2451599.77 ±0.08 2451600.03 ±0.08 

62.2 59.0 ± 0.6 60.8 ± 0.6 



For each of the above three datasets, the average and the stan- 
dard deviation were calculated. It is clear from the results, 
which are listed in Table[5] that the secondary's RVs are in 
much better agreement than those of the primary. The larger 



{Av(t)} t = {Av, 



(t)} t U {Av secondary (*)}*• 



average difference of the {At 



primary 



(t)}t dataset is due to the 



primary's line profile being so strongly rotationally broadened. 
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Table 3. Statistical properties of the differences between the 
final RV values derived with iteration scheme II and I. See text 
for more information on the three datasets. 



Dataset 


Average (kms 


x ) Standard deviation (km s L ) 


{Al) primary (t) }t 


2.37 


1.02 


{At? secondary (^) }t 


0.26 


0.27 


{A«(t)} t 


1.32 


0.75 



This implies that the orbital solution derived by means of only 
the secondary's final RVs is very stable regardless of the ap- 
plied iteration scheme, while the orbital solution derived on the 
basis of only the primary's final RVs is somewhat less stable 
and, hence, causes the small difference between the final results 
of iteration scheme I and II. We consider the results obtained 
with iteration scheme II more reliable because they are in better 
agreement with the interferometrically derived orbital solution 
on the one hand and, although it is not clear in the rounded val- 
ues listed in Table|2j the uncertainties in the derived parameters 
derived with scheme II are all smaller. 

We also note that both iteration schemes lead to smaller er- 
ror bars on the derived orbital parameter values than the ones 
obtained in Paper I. FiguresQ and [8] show that the uncertain- 
ties in the system velocity, the semi-amplitude of the primary's 
orbit, and the longitude of periastron have been significantly 
lowered by the iterative process. The application of iteration 
scheme II maintains or even improves the compatibility with 
the interferometric results (see Fig.|8|l. Therefore, we conclude 
that the iterative process results in a significant improvement of 
the orbital solution. 

The final primary's (open dots) and secondary's (filled dots) 
RV values that were obtained with iteration scheme II are 
shown in Fig.[9] The best combined fit to these datasets is 
denoted as a full (primary's orbit) or dashed (secondary's or- 
bit) line. The fit is satisfying compared to the one obtained in 
Paper I. 

Finally, we applied KOREL again with this new, iteratively 
derived orbital solution. Indeed, in theory this can result in 
improved versions of the disentangled profiles of both com- 
ponents and, therefore, allow an iteration process on a higher 
level. However, the best disentangled profiles that came out of 
the KOREL analysis were hardly distinguishable from the ones 
shown in Fig.|4]and, hence, there is no point in repeating the it- 
eration process with these "new" disentangled profiles. In par- 
ticular, the small bump in the centre of the primary's disentan- 
gled profile in Fig.|4]did not disappear. It is due to the imperfect 
averaging over the oscillations of the primary by KOREL. This 
is not surprising in view of its complex multiperiodic high- 
degree oscillations, which we cannot unravel perfectly from 
our data (see Sect. 4). 

3. Physical parameters of the components 

The physical parameters of /3Cen derived in PapeiTI were ob- 
tained without taking into account the systematic effects de- 
scribed here. Their values and errors must clearly be revised. 
The error estimates were optimistic as they were derived from 
systematically underestimated RV values. 



We first checked that the refined value of the orbital period 
(which was fixed in deriving the interferometric orbital param- 
eters in Paper II) does not change the values of the orbital in- 
clination and angular semi-major axis. We then followed the 
same strategy for the computation of the individual compo- 
nent masses as in Paper II, using the refined values of the or- 
bital period and the semi-amplitudes of both components from 
scheme II, the orbital inclination from the interferometry given 
in Paper II, and the mean of the spectroscopic and interfero- 
metric values for the eccentricity. We also need to take into 
account the standard errors of these quantities. This is not so 
straightforward for the parameters resulting from our iteration 
schemes because their la error from FOTEL listed in Table[2]is 
necessarily an underestimation of the true error. Indeed, these 
errors were derived under the assumption that the disentan- 
gled profiles resulting from KOREL are error-free, which is not 
the case. We are, unfortunately, unable to propagate the un- 
certainty induced by the disentangling properly because the 
current version of KOREL does not provide us with error es- 
timates. Moreover, our schemes implicitly assume that KOREL 
appropriately treats the effects of the oscillations as random 
noise in computing the disentangled profiles. For this reason, 
we adopt a very conservative approach for the error propaga- 
tion and use 2a errors rather than those listed in Table|2]in the 
derivation of the physical parameters of the system. Following 
the approach of Paper II, this leads to Mi = 11.2 ± 0.7 M Q 
and M-2 = 9.8 ± 0.7 M . Finally, this explains our choice for 
the revised nomenclature of the primary and secondary. 

We attempted to refine these estimates by using all the addi- 
tional observational information we have at our disposal, such 
as the CORALIE echelle spectra. We first estimated the effec- 
tive temperature and gravity of the two components by follow- 
ing the procedure outlined in Uytterhoeven et al. (2005) for 
the double-lined binary k Sco, i.e., by merging theoretical line 
profiles of H, He, and Si lines with the appropriate flux ratio 
according to NLTE predictions made from the latest version 
of the FASTWIND code (Puis et al. 2005), after using the or- 
bital RVs to shift the profiles. This led us to the conclusion that 
both components have T off = 24 000 ± 1 000K and log g = 
3.4 ± 0.3. The large uncertainty in the gravity stems from the 
difficulty in achieving a proper normalisation of the spectra 
near the Balmer lines. Since we find the two components to 
have equal T e ff and log g within the uncertainties, it is possible 
to compute photometric estimates of these quantities from mul- 
ticolour photometry. We did this from Geneva measurements 
of /3Cen at our disposal and find T eS = 26 500 ± 500K and 
log g — 3.7 ± 0.2 assuming equal components. This leads 
us to a safe broad range of T e s = 25 000 ± 2 000K and 
log g = 3.5 ± 0.4 for both components. 

We subsequently scanned the very extensive database of 
main-sequence stellar models published by Ausseloos et al. 
(2004), which have a range in mass from 7 to 13 M Q in steps of 
0.1 M Q and a range in Z from 0.012 to 0.030 in steps of 0.002, 
for each of the three values of the core overshooting parameter 
of 0.0, 0.1, and 0.2 expressed in local pressure scale heights. 
The models have X — 0.70 and the solar mixture of Grevesse 
et al. (1996). For a description of the input physics, we refer to 
Ausseloos et al. (2004). We scanned this database requesting 
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Fig. 9. The orbital curves obtained by means of iteration scheme II. The open and filled dots represent the RVs of, respectively, 
the primary and secondary components around the centre of mass. The full and dashed lines represent the best-fitting orbits of 
the primary and secondary, respectively, according to the parameters listed in the column labelled "Scheme II" in Table[2] 



that the masses, effective temperatures, and gravities of /3 Cen's 
components lie in the appropriate ranges and that the age of 
the components must be equal to within 1%. This leads us to 
acceptable ranges for the masses of Mi G [10.6, 10.8] M Q and 
M 2 S [10.2, 10.4] M®, and an age e [13.5, 14.7] million years. 

As an a posteriori check, we computed the allowed mass 
ratio resulting from the brightness ratio obtained from the 
interferometry in PapeiTI and the mass-luminosity relation 
log(L 2 /Li) = (3.51 ± 0.14)log(M 2 /Mi) (Griffiths et al. 
1988). We thus find the condition M 2 /Mi e [0.95,0.97], 
which is fulfilled by our solutions for the masses resulting from 
the spectroscopy and interferometry refined by the modelling. 

Finally, the systematic errors in the semi-amplitudes also 
call for a re-evaluation of the dynamical parallax given in 
PapeiTI. Following the same approach as in Paper II and using 
2cr errors for the spectroscopic elements, we find ir = 9.3±0.3, 
resulting in a distance of 108±4 pc. 

4. Analysis of the line-profile variability 

Challenging aspects of massive star asteroseismology are the 
detection of numerous frequencies and their mode identifica- 
tion on the one hand, and the derivation of the fundamental pa- 
rameters of the targets (T c g, log g, M) with high precision on 



the other hand. We succeeded in the latter, and make an attempt 
to tackle the former challenge now. 

4. 1. Frequency analysis 

Our aim is to find the timescales associated with the short-term 
line-profile variability and to connect each timescale with the 
component to which it belongs. This is by no means straight- 
forward because the short-term variations have a significantly 
lower amplitude than the orbital variations. Moreover, the line 
profiles of both components are fully blended with each other 
at all orbital phases. This required a specific non-standard anal- 
ysis, the details of which are available in Ausseloos (2005). 
Here, we present only a concise summary of the results. In par- 
ticular, we point out that the complexity of the profile varia- 
tions due to the presence of moving subfeatures (see Ausseloos 
2005, p. 28, Fig. 2.1 for examples and Fig.|5|l does not favour a 
standard radial-velocity analysis, but requires a search for fre- 
quencies across the whole width of the profiles. 

We perform an analysis of the intrinsic variability of the 
primary by using a two-dimensional (2D) frequency analysis 
method first introduced by Gies & Kullavanijaya (1988) and 
later defined as the Intensity Period Search (IPS) by Telting et 
al. (1997). We used a 2D version of the Lomb-Scargle method 
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Fig. 10. Grey-scale representation of the power spectrum ob- 
tained at each wavelength position (given in units of the vari- 
ance) across the Silll A 4552.6 A line profile by the application 
of the 2D Lomb-Scargle method to all CES spectra. 



(Scargle 1982) for the time series of normalised flux values 
at each wavelength across the profile. As we found a value of 
0.0045 c d _1 for the half-width at half-maximum of the central 
peak of the window function based on all CES data, and as the 
2D frequency analysis is rather time consuming, we adopted 
a frequency step of 0.001 cd -1 in a first stage. After having 
identified the main peak and its aliases, we recomputed the pe- 
riodograms with a factor 10 smaller frequency step around the 
dominant peak and its aliases to check if the results remained 
valid, which was always the case for the relevant frequencies 
mentioned below. 

There exists a small wavelength range centred around 
4549 A, at which only the primary's line profile is present. We 
first applied the 2D Lomb-Scargle method to the dataset com- 
prising all the CES spectra. A graphical representation of the 
results is shown in Fig.^| This figure reveals power excess 
in the wavelength range between 4550 A and 4557 A, which 
is visible at frequencies below 3 cd _1 . This is caused by the 
shift of the secondary's line profile due to its orbital motion. 
However, clearly visible peaks occur between 5 and 9 c d _1 in 
the power spectra at wavelengths between 4549 A and 4550 A. 
This excess power can only be due to the primary, proving that 
this component has short-term periodic variability. 

To find the frequencies of the short-term variability of 
the primary with better significance, we removed the higher- 
amplitude variability due to its orbital motion around the cen- 
tre of mass. For this, we used the secondary's disentangled 
line profile in combination with the orbit. We then computed 
the periodograms at each wavelength to construct a 2D peri- 
odogram. The results for the CES spectra are shown in Fig.ll II 
(we omitted the CORALIE spectra for this plot due to their 
larger noise level). We clearly reveal excess power across the 
whole line profile of the primary at a frequency near 6.4 c d _1 
and its aliases. 

We summed all the ID periodograms (for both the CES and 
CORALIE spectra) over the range 4549 and 4556.5 A (Fig. [HI. 
The drawback is that this analysis method removes the sec- 
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Fig.ll. Same as Fig. [TO] but after subtracting the secondary's 
disentangled line profile and correcting for the primary's orbital 
motion. 
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Fig. 12. Power spectra resulting from the addition of all 
the ID periodograms over the wavelength range between 
4549 and 4556.5 A. The panels correspond to different stages 
of prewhitening: original periodogram (upper panel), af- 
ter prewhitening with 7.415 cd -1 (middle panel), and after 
prewhitening with 7.415 and 4.542 c d _1 (lower panel). 



ondary's long-term variability, but not its short-term variability. 
Therefore, the sum of all ID periodograms shown in Fig. 1 121 
contains a mix of peaks due to the primary's and secondary's 
short-term variability, if any. The upper panel unquestionably 
reveals 7.415 c d _1 as the dominant frequency, although rather 
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strong aliasing still occurs for the whole dataset. We have 
shown above that this frequency belongs to the primary. The 
middle panel of Fig.lT2"lshows the periodogram that is obtained 
after prewhitening the original data with 7.415 cd _1 . It sug- 
gests 4.542 c d _1 , or one of its aliases, as the second frequency. 
After prewhitening with 7.415 and 4.542 c d~\ the power spec- 
trum (lower panel) is dominated by peaks at lower frequencies. 
The highest peak in the interval [3,10] cd -1 appropriate for 
/3Cephei stars occurs at 4.407cd~ 1 , but it is not clear at this 
stage whether this is another intrinsic frequency. We conclude 
that there is evidence for unexplained additional power. 

A study of the secondary's intrinsic temporal behaviour is 
by no means straightforward. Due to their large width, the pri- 
mary's line profiles extend over a wavelength range that com- 
pletely includes the range spanned by the secondary's line pro- 
files during nearly all orbital phases. This makes it very difficult 
to unravel the secondary's line-profile variations, if any, from 
the ones of the primary. We applied the following procedure: 
for each original spectrum, the primary's Silll A 4552.6 A dis- 
entangled line profile was shifted according to the correspond- 
ing orbital velocity and subsequently subtracted from the orig- 
inal spectrum. We then computed a 2D Scargle periodogram 
as explained above. We subsequently added the power across 
subintervals of the total wavelength range [4547, 4557] A. At 
the same time, we determined the extent of the secondary's 
line within that wavelength range (see Ausseloos 2005, p. 72, 
Fig. 2.32 for more details). If the power is only significant in 
the wavelength range spanned by the secondary, this is consid- 
ered as a strong indication that the corresponding frequency is 
associated with that star. 

We carried out several tests by considering different sub- 
datasets whose power distribution for candidate frequencies 
was calculated over different subintervals in wavelength. From 
all these tests, we conclude that also belongs to the primary 
and that no evidence of short-term variability in the line profiles 
of the secondary was found in our dataset. It has to be stressed 
that all previously published period analyses of (3 Cen (Breger 
1967; Shobbrook & Robertson 1968; Lomb 1975; Robertson 
et al. 1999; Ausseloos et al. 2002) can no longer be trusted as 
they all neglect the pulsations of the primary and assume that 
the component with the deeper and narrower line profiles un- 
dergoes the short-term variability. 



4.2. Mode identification 

An attempt was made to identify the modes of the primary. 
Given its strongly rotationally-broadened profiles, the Doppler 
Imaging method was used. This method was introduced by 
Gies & Kullavanijaya (1988) for the B0.7III star e Per, but sev- 
eral authors have elaborated on it since (see, e.g., Telting & 
Schrijvers 1997 and references therein). Telting & Schrijvers 
(1997) took a major step forward by making a large Monte- 
Carlo simulation study from which they derived linear relation- 
ships between the degree £ and the blue-to-red phase difference 
f of an observed frequency / on the one hand, and between 
the azimuthal number m and the phase difference of the first 
harmonic A^>2f on the other hand. The errors of the estimates 
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Fig. 13. Phase behaviour across the profile calculated on 
the basis of all CES and CORALIE data for (/i,/ 2 ) = 
(7.415, 4.542) cd _1 . The upper (middle) panel shows the phase 
behaviour of /i (/2). The average line profile is given in the 
lower panel. 



of £ and m are, respectively, one and two. Telting & Schrijvers 
(1997) also verified that the method can handle multiperiodic 
line-profile variations. 

We applied this method to different combinations of alias 
frequencies (/i,/a) with f\ = 7.415 c d _1 or one of its adjacent 
aliases, and f 2 = 5.546 or one of its adjacent aliases and the 
results are the same for each of them. The resulting phase be- 
haviour for (/i,/ 2 ) = (7.415, 4.542) cd" 1 is shown in Fig.[l3] 
We obtain smooth phase distributions across the line profile 
which allow a reliable application of Telting & Schrijvers' 
(1997) linear relations to estimate I. For the first frequency 
fi = 7.415 cd~\ we can read off a blue-to-red phase dif- 
ference f 1 of 57r radians, irrespective of the value of fa. 
This implies that £\ G [4, 7]. The phase diagram of 4.542 cd -1 
leads to a phase difference A 1 ^ f 2 e [3.5, 3.6] 7r radians, hence 
£2 G [3,5]. These results explain why variations in ground- 
based photometry are absent because £ values above two lead 
to strong cancellation across the visible stellar disk in such data 
(see, e.g., Dziembowski 1977). 

We scanned the database of seismic (3 Cephei star models 
and their oscillation frequencies computed by Ausseloos et al. 
(2004) once more, considering the tight limitations on the phys- 
ical parameters of (3 Cen's primary. Despite the narrow allowed 
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range in the mass, the effective temperature, and the age of the 
primary, we could find numerous predicted modes with a fre- 
quency near f\ or f 2 for the allowed models. As we have no 
definitive mode degree for the two frequencies, nor any esti- 
mate of their azimuthal orders, we cannot refine the physical 
parameters of the primary from the oscillations at this stage. 

5. Conclusions 

We have shown a systematic error to occur in the semi- 
amplitudes of the velocity curves, due to an underestima- 
tion of radial-velocity values in SB2 spectroscopic binaries 
with merged spectral line profiles. We provide methodological 
schemes to solve for this systematic error. They are based on 
spectral disentangling by means of the KOREL code (Hadrava 
1997). We suggest that these schemes be used in any future 
analyses of SB2s whenever their profiles can be successfully 
disentangled. 

In the case of (3 Cen, the systematic underestimation of the 
spectroscopic orbital semi-amplitudes led to an underestima- 
tion of the component masses of about 10%. We refined the 
component masses of this massive binary by application of our 
analysis schemes to the available high-resolution spectroscopy, 
and by combining the spectroscopic results with interferomet- 
ric measurements across the orbit, leading to a precision of 6%. 
The accuracy was further improved by stellar modelling taking 
into account an extensive database of stellar evolution models 
with wide ranging values of the mass, Z, and core convective 
overshooting. In this way, we find the component masses of 
/3Cen to be M x = 10.7 ±0.1 Mq and M 2 = 10.3 ±0.1 Mq 
and its age to be 14.1 ± 0.6 million years. These mass esti- 
mates turn out to be fully compatible with the mass-luminosity 
relation. The fact that we find (3 Cen to have passed less than 
half of its main-sequence lifetime is compatible with its high 
eccentricity and suggests that both components were formed 
together, rather than having undergone a tidal capture. The ab- 
sence of an IR excess (Aerts et al. 1999), and of Ha emission in 
the CORALIE spectra, exclude the stars still being in their pre- 
main-sequence phase. The a posteriori conclusion that the de- 
rived fundamental parameters of the components of the system 
fulfill the tight mass-luminosity relation provides confidence in 
our high-precision estimates of the masses. The determination 
of the distance to (3 Cen is also affected by the previous system- 
atic underestimation of the radial-velocity values and has been 
re-determined from its dynamic parallax to be 108±4pc. 

Next, we performed an in-depth line-profile analysis. All 
previous frequency analyses were focused on the line-profile 
variations of the component producing the deeper lines (sec- 
ondary) and claimed frequencies for this component. We have 
given compelling evidence that it is actually the component 
producing the broader lines (primary) which is pulsating, while 
we did not find any periodic variability that we could link to the 
secondary. As our dataset is by far the most extensive one so far 
used for a spectroscopic analysis of (3 Cen, the pulsational na- 
ture of the secondary should be regarded as unproven. If the 
secondary oscillates, its amplitudes must be much smaller than 
those of the primary because the variations of the combined line 



profiles of both components are dominated by the variations of 
the primary's line profile. 

We detected two frequencies in the primary's line-profile 
variations by means of the 2D Lomb-Scargle method, but we 
were not able to fix their value unambiguously due to aliasing. 
Notwithstanding the aliasing effect, we were able to restrict the 
value of the degree I to [4, 7] and [3, 5] for the first and second 
mode, respectively. The detection of only two frequencies and 
the lack of unique mode identification prevented an in-depth 
seismic study of the star, despite the fact that its fundamental 
parameters are so tightly constrained by the binarity. The only 
conclusion we could draw in this respect is that standard stellar 
models predict frequencies that are fully compatible with the 
two detected ones. 

Very few accurate masses of (3 Cephei stars are available, 
notable exceptions being those with a seismic mass estimate 
(HD 129929: Aerts et al. 2003; 16 Lac: Thoul et al. 2003; v Eri: 
Pamyatnykh et al. 2004 and Ausseloos et al. 2004). The masses 
we derived here for (3 Cen, together with the estimates for its ef- 
fective temperature, gravity, and age, constitute a fruitful start- 
ing point for future seismic analyses of this massive binary. A 
distinct short-event photometric variation of (3 Cen with am- 
plitude of about 0.04 mag was observed by Balona (1977) on 
one night. No period could be derived in these data, how- 
ever. The lack of other claims of photometric variability of the 
brightest among all (3 Cephei stars, despite observational efforts 
(L.A. Balona, private communication), is nicely explained by 
our detection of pulsation modes with a high degree (£ > 3). 
This clearly points out that one cannot hope to find a com- 
plete frequency spectrum of the p-modes in (3 Cephei stars from 
ground-based photometry alone. The same conclusion was re- 
cently drawn for the fast rotator ( Oph on the basis of high- 
resolution spectroscopy and MOST space photometry (Walker 
et al. 2005). It is clear that high-precision photometric data 
from space are necessary to achieve a seismic interpretation of 
/3Cen. 
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